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ABSTRACT 

Recent discoveries of strongly misaligned transiting exoplanets pose a challenge to 
the established planet formation theory which assumes planetary systems to form and 
evolve in isolation. However, the fact that the majority of stars actually do form in star 
clusters raises the question how isolated forming planetary systems really are. Besides 
radiative and tidal forces the presence of dense gas aggregates in star-forming regions 
are potential sources for perturbations to protoplanetary discs or systems. Here we 
show that subsequent capture of gas from large extended accretion envelopes onto 
a passing star with a typical circumstellar disc can tilt the disc plane to retrograde 
orientation, naturally explaining the formation of strongly inclined planetary systems. 
Furthermore, the inner disc regions may become denser, and thus more prone to speedy 
coagulation and planet formation. Pre-existing planetary systems are compressed by 
gas inflows leading to a natural occurrence of close-in misaligned hot Jupiters and 
short-period eccentric planets. The likelihood of such events mainly depends on the 
gas content of the cluster and is thus expected to be highest in the youngest star 
clusters. 

Key words: hydrodynamics — planets and satellites: formation — planet-disc in- 
teractions — planet-star interactions — protoplanetary disks — open clusters and 
associations: general 



1 INTRODUCTION 

The discoveries of retrograde or strongly m isaligned tran- 
siting exoplanets l i ke WASP- 17b and others (JHebrard et all 



l2008V'Gillon"2009'; 'Johnson et al."2009'; Narita et al. 2009; 
Pont et al. 2009, 2010; V^ inn et al. 2009a,b; Anderson et al. 
I2OIOI : iTriaud et al.l I2OIOI I pose a challenge for established 
planet formation theories. These classically assume the birth 
of planetary systems out of collapsing protostellar cloud 
fragments. While contracting from sub-parsec to millipar- 
sec (or hundreds of AU) scale the conservation of its initial 
angular momentum causes the cloud to spin and flatten, 
while part of its mass settles towards the center, eventually 
igniting the hydrogen burning in the stellar core. Given such 
a protostar with a circumstellar disc that contains a certain 
amount of dust, solid particles collide and stick together, 
forming larger particles. This coagulation leads to bodies 
(cores) of sufficient mass to accrete surrounding material 
(dust and even gas) by gravity until no more disc mate- 
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rial is available. This mechanism is probably the dominant 
process of the formation of the Solar System, including the 
Earth (Morda sini et al." 2008). As an alternative, gravita- 
tional instability as a forming mechani sm for gian t planets 
and brown dwarfs is being discussed (jBossI 1 20041 ) , and at 
least for brown dwarfs it has actually been shown to work 
( Stam atellos et al.1 [JOOTal : iThies et al.ll2010l ). More recent 
variations of these models assume gravitational instability 
of the solid phase only or dust trapping in vortices as a 
speed-up for coagulation. 

Being isolated from any external perturbation, every- 
thing inside this protostar-disc system spins in the same 
direction, thus the forming planets orbit their host star the 
same way. This simple model has been severely questioned 
by the findi ngs mentioned above. The transiting exoplanet 
WASP-17b (lAnderson et al.ll2010l ). for example, has been 
found to have a sky-projected inclination of the orbital plane 
normal against the stellar spin axis of I48.5I42 degrees. A 
number of other mis-aligned transiting exoplanets have been 
discovered JHebrard et al."2008': 'Gillon"2009'; 'John son etllT 
20091: [Narita et al.ll2 009: Pont et al. 2009, 2010; Wi nn et al 



2009al lbl: lAnderson et al.l |2010|) • These discoveries suggest 
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that misalignment between planetary orbits and the spin 
of their host star (hereafter called spin-orbit misalignment) 
are quite common, at least among close-in transiting planets 
for which the spin-orbit alignme nt can be measured by using 
the R ossiter-McLaughlin Effect (|Rossitei]ll923 : iMcLaughlinI 
I1924I I. The nearby v And planeta ry system exhibits ev en 
mutually inclined planetary orbits (jMcArthur et al.ll2010h . 
Such planetary orbits cannot be explained by the stan- 
dard planet formation model according to which the con- 
servation of angular momentum forces any coplanar system 
to remain coplanar forever. There have been recent sugges- 
tions that may shed some ligh t into the pos s ible m echanism. 
iFabrvckv fc Tremaind (|2007l ): iNaoz et al.1 (|201in deduce a 
formation scenario for misaligned "hot jupiters" through 
long-term mutual perturbations of two planets (or an inner 
planet and an outer brown dwarf; for simplicity, we use the 
term "planet" for all substellar companions in this paper) 
orbiting the same star. Given an initial mutual inclination 
between 40° and 140° the Kozai-Lidov mechanism leads to 
secular oscillations of eccentricity vs. inclination. Tidal fric- 
tion circularises the orbit of the inner planet if its periastron 
falls below a few stellar radii, eventually leaving the planet 
on a close-in orbit with near-random inclination with re- 
spect to the stellar spin. This scenario, however, requires a 
sufficient initial mutual inclination of both planets. 

Another possible mechanism is tilting the spin axis of 
the star. Direct stellar encounters within a stellar radius 
(about 0.01 AU) can be effectively ruled out, b ecause they 
are improbable, even within stellar clusters (JThies et al.l 
I2OO5I : iThies fc Kroupal l2007l ) and would destroy the disc. 
Secular transfer of angular momentum between the proto- 
stellar spin and the protoplanetary disc v ia magnetic fields 
have been investigated by Lai et al.l (|201lh . 

Multi-stage or episodic accretion of circumstellar mate- 
rial may provide another viable mechanism for misaligned 
planetary systems. The underlying idea is that accretion 
from different sources (i.e. gas filaments, accretion en- 
velopes) from different directions may sever the classically 
assumed spin-orbit correlation. The star itself may get its 
angular momentum from a different accretion event than the 
bulk of circumstellar material. Bate et al. (2010) have anal- 
ysed t urbulent accretion events using data from iBate et al.l 
(|2003l ). They deduce a high frequency of misaligned plan- 
ets as a consequence of multi-directional accretion. How- 
ever, due to the limited simulation data available, they 
were forced to make simplified assumptions to what hap- 
pens within the close vicinity of a star. Our work aims to 
take a deeper insight into such scenarios by self-consistent 
calculations of close interactions of a circumstellar disc with 
encountered new material, and by also treating the effects 
on pre-existing planets. 

As IPfalzner etHI (|2008l l: lOlczak et all (|2010h have 
shown, encounters between stars and pre-/protostellar ob- 
jects do occur frequently in dense stellar environments and 
may lead to significant acc retion bursts. Althou gh the pre- 
cise mechanisms proposed (|Pfalzner et al.ll2008l ) differ from 
ours, their calculations emphasise the importance of en- 
counter events to accretion processes. 

In this paper the effects of gas capture onto a pre- 
existing circumstellar disc are studied. In Section [2] we depict 
the physical scenario while Sect ion [3] briefly describes the nu- 



merical methods used in our computations. The results are 
presented in Section |4l 



2 MODEL 

In our model, both the initial mutual inclination as well as 
the inclination with respect to the stellar equator can be 
naturally explained by multi-stage accretion during the for- 
mation process of the planetary system. The scenario can be 
described as follows: First, a star with a protoplanetary disc 
(PPD) is forming the classical way, i.e. out of a collapsing 
spinning cloud fragment. Soon after this the protostar passes 
another region in the star cluster that contains dense gas, 
e.g. the protostellar cloud fragment or accretion envelope of 
another star, hereafter called the target. Accordingly, the 
star which encounters the target and its disc are referred to 
as the bullet star and the bullet disc, respectively. In our cal- 
culations, a massive circumstellar disc around another star 
is used as the target. While grazing the outer regions of this 
target the star captures material from it at an angle that 
depends on the encounter orbit and the arbitrarily chosen 
orientation of the target and the bullet disc. Since there is no 
correlation between the orientations of the bullet disc, the 
target and the encounter orbit, this configuration reflects an 
uncorrelated or random in-flow of material. The subsequent 
evolution mainly depends on the amount and infall angle 
of this additional material, and also on the protoplanetary 
evolution stage during the capture. 

In our calculations we focus on the case of an encounter 
before any planet has formed in the passing star's disc. It 
is thus a purely gaseous interaction computed via smoothed 
particle hydrodynamics (SPH). In our model, the star is a 
sun-type star with a low-mass PPD (Md < 0.1 Mq, radius 
< 100 AU), passing through the target disc (0.5 Mq). 

The initial conditions for the P PD are taken from 
IStamatellos fc WhitworthI (|2008l . l2009l l . The disc model has 
power-la w profiles for temperature , T, an d surface density, 
S, from Stamatellos & WhitworthI (|2008), both as a func- 
tion of the distance, R, from the central star 
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where 1.5 < ^s < 1.75 and 0.5 < qr < 1.0. Too = 10 K is a 
background temperature to account for background radia- 
tion from other stars within the host cluster. The radiation 
temperature of the stars are estimated assuming a mass-to- 
luminosity relation of L oc M^ for low-mass stars. 

The disc is initially populated between 4 and 40 AU. 
Preceding the actual encounter the disc is allowed to "settle" 
during the approach time of about 5000 years (equivalent to 
20 orbits at the outer rim) , erasing artefacts from the initial 
distribution function. During this settling the disc smears 
out at the inner and the outer border extending its radial 
range from a few AU up to > 60 AU. The actual surface 
density after the settling and thus the mass content within 
a given radius is therefore slightly smaller than in the initial 
setup. 
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Similarly, the target is set up as a massive extended disc 
with mass of 0.48 Mq and a radius of ~ 400 AU, as well as 
a 0.75 Mq star to keep the gas gravitationally boun d. The 
setup is therefore equivalent to those used in Thie s et al.l 
([2OIOI), where triggered fragmentation inside this envelope 
upon tidal perturbations have been studied, and is indeed 
realistic as long as encounters with large circumstellar discs 
are considered. However, future models will also use different 
setups (protostellar cores, dense gas filaments etc.) as the 
target. 

Both systems are initially set on a hyperbolic orbit 
around their centre of mass with the discs being inclined 
mutually as well as with respect to the orbital plane. The 
periapsis is set to 500 AU initially. The eccentricity is set 
to 1.1, corresponding to a pre-encounter relative velocity of 
about 0.7 km~^. The periapsis relative velocity is 2.9 km~^. 
The hyperbolic orbit is inclined 15 degrees against the tar- 
get disc plane using Eulerian z-x-z rotation angles 0, 15°, 0, 
while the plane of the bullet disc (initially in the x-y plane) 
is tilted by 0, 135°, 60°. 

The eccentricity value is in agreement with th e likeli- 
hood estimation of encounter parameters by Olcz ak et al.l 
(|20ld V In addition, encounters at distances of the order of a 
few hundred AU are likely to happen in clusters comparable 
to th e Orion Nebula Cluster JThies et al. II 20051 : 1 Olczak et al.1 
I2OIOI), thus our scenario is rather likely than exceptional. 
The actual periapsis and eccentricity may differ slightly due 
to tidal angular momentum transfer, but these effects on the 
orbit are negligible at that time. After the passage, when 
both stars have reached a sufficient distance {^ 1000 AU) 
again, the gas particles around the bullet star (i.e. the re- 
mains of its own disc and any captured gas) are separated 
from the system and analysed in the centre-of-mass system 
of the bullet star. 



3 NUMERICAL METHOD 

All computations were performed by using the well-tested 
sm oothed-particle h ydrod ynamics (SPH) code dragon 
by I Goodwin et al.l (|2004|) including the radiati ve heat 
transfer extension by IStamatellos et all (|2007bl l. Most 
of th e numerical parameter s et tings have been adopted 
from I Stamatello s et al.l (|2007al l: IStamatellos fc Whitworth 
((20091). The radiative-transfer algorithm is not strictly re- 
quired in our model, but due to its high numerical efficiency 
it does not significantly slow down the calculations, so we 
left it in to keep i t consistent with previous calculations by 
I Thies et al.l (|2010l V All gas particles have the same mass 
and equation of state. The artificial viscosity parameters are 
a = 0.1 and /3 = 0.2, i.e. lower then th e usually used values 
(1 and 2, respectively, Monaghan 199Z), to keep the disc dis- 
persion low. Stars and planets are treated by sink particles 
which are set in the initial conditions, but may also form if 
the local volume density exceeds 10~^ g cm~^ (not expected 
in the current model, however). The sink radius is chosen 
as 0.5 AU while the sink masses are the masses of the stars 
involved, i.e. 1.0 Mq for the bullet star and 0.75 Mq for the 
host star of the target disc. Any gas particle that becomes 
gravitationally bound to a sink within less than the sink 
radius is treated as being accreted by the sink. The code 
currently does not treat magnetic field nor dust. 
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Figure 1. Snapshot of a bullet star with a 0.1 Mq circumstellar 
disc passing through a massive extended (0.5 Mq, 400 AU) target 
disc on an inclined orbit. While the original disc is partially trun- 
cated, partially condensed to the center, and somewhat tilted, 
the bullet star accretes additional material from the target disc, 
forming an inclined annulus of gas around it (frames b to d). The 
initially strongly inclined structures progressively align each other 
within 30 000 years after the encounter (frames e, f). The time 
stamp in each frame refers to the time of the encounter. 



The numerica l model used here is the same as used in 
iThies et al.1 (|2010h . using 250 000 SPH particles for the tar- 
get while the PPD of the passing star has been set up by 
25 000 to 50 000 particles, for disc masses between 0.05 to 
0.1 Mq, respectively. 



4 RESULTS 

4.1 General effects 

Figure [T] shows six snapshots of the encounter and subse- 
quent gas capture. The view is centred on the passing star 
(the "bullet star") at face-on direction towards the unper- 
turbed disc. The time stamp refers to the moment of the pe- 
riastron passage. Frame (a) shows the unperturbed disc. Its 
core has a radius of about 40 AU while the low-density outer 
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Figure 2. The angular momentum unit vectors of the particles of 
the material around the bullet star in spherical coordinate angles, 
1? (latitude) and (^ (longitude), i.e. projected on the unit sphere. 
The reference frame is adjusted so that the normal vector of the 
bullet disc points at (<^,^) = (0, 0). The topmost two frames show 
the angular momenta initially and at the time of the closest ap- 
proach. While the original bullet star disc is altered only slightly 
there is clearly another population of particles settling in an annu- 
lus which is tilted against the initial disc plane. After formation, 
both structures progressively align to each other, and the angular 
difference between both degrades (subsequent frames). The total 
angular momentum vector of each disc is marked by a cross, while 
the open circle indicates the initial angular momentum vector of 
the bullet star disc before the encounter (being in fact nearly 
identical to that at the moment of the encounter, t = 0). This 
also denotes the direction of the bullet's star spin vector. The 
time stamp in each frame refers to the time of the encounter. 



regions extent to about 100 AU. Shortly after the encounter, 
as shown in frame (b-d), the captured material starts to 
form an annulus around the star while, at the same time, 
the original bullet disc is partly truncated, with the inner 
region being condensed. The persistent quadrupole forces 
continue to alter the disc's orientation and eventually tilt 
it into edge-on orientation (frames e and f). At the same 
time, the orientation of the annulus is only slightly changed 
due to its larger radius and thus higher moment of inertia. 
As clearly visible, both the captured annulus and the disc 
are progressively aligning each other, eventually forming a 
combined disc. The dynamical evolution of the mutual orien- 
tation and the total tilt of the original bullet disc is depicted 
in Figures [2] and O Here, the SPH particles have been sepa- 
rated into groups sharing a similar angular momentum unit 
vector and thus belonging to the same rotational structure. 
Particles which deviate more than 15 degrees from the bulk 
angular momentum as well as unbound or loosely bound 
particles beyond 100 AU from the bullet star are neglected. 
The mutual inclination of both structures is declining from 
initially 130 degrees to near zero degrees within 35 kyr. On 
the other hand, the original disc is tilted by nearly 100 de- 
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Figure 3. The mutual angular tilt between the pre-existing disc 
of the bullet star and the annulus of captured gas, as well as the 
net angular tilt of the pre-existing disc as a function of time. While 
being almost retrograde wrt. each other before the encounter, 
the two structures progressively align with each other over time. 
25 kyr after the encounter the mutual tilt has degraded to about 
35 degrees. At the same time, the pre-existing disc has been tilted 
wrt. its initial orientation by about 90 degrees. It has become 
retrograde wrt. the bullet star's spin after about 27 000 yr. 



grees, i.e., if being initially aligned to the stellar equator it 
is now retrograde. 

During the encounter, typically around 10 to 30 Mj of 
gas is captured by the bullet star (Figure [T] frames b-d; 
Section 132)) • In the cases where a protoplanetary disc (here- 
after called the original disc) is already present, two general 
effects can be distinguished: (i) During the capture process, 
the original disc is compressed to about one-half or one-third 
of the initial radius, accordingly increasing the surface den- 
sity, (ii) The captured material settles in an outer annulus 
the inclination of which depends on the orbital orientation 
and the orientation of the target disc relative to the en- 
counter orbit. 

The orientation of this annulus, shown in Figure [2] in 
spherical coordinates, is not constant, but changes over time. 
Within about a few thousand years after formation the an- 
nulus and the bullet disc have considerably lowered their 
mutual inclination (third and fourth frame), as shown in 
Figure [3] A few ten thousand years after the encounter, 
they have almost aligned to each other, eventually form- 
ing a single circumstellar disc. Consequently, the orientation 
of this disc differs considerably from its initial orientation 
and therefore from the stellar equatorial plane. It should be 
noted that this tilt may be muc h lower if a steady and h o- 
mogeneous gas flow is assumed (|Moeckel fc Throodl2009[ ). 



4.2 Captured mass 

The amount of captured mass is an important quantity to 
judge its impact on planet formation. In our calculations, 
typically between 0.01 and 0.02 Mq is captured from the 
target disc, i.e. about one third of the initial bullet star disc 
mass. Figure |4] depicts the time evolution of the masses of 
the bullet star disc and the captured annulus. As shown here, 
the capture of material from the target disc begins quickly 
about the closest encounter, and is almost finished within 
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Figure 4. The mass of the pre-existing bullet star disc (solid 
curve) and of the captured annulus (dashed curve; see Section 
14. 1|) as a function of time. The peaks in the solid curve at about 
the moment of the closest encounter are due to the perturbations 
of the pre-existing disc caused by the impact of in-flowing gas and 
tidal forces by the target mass. The capture process starts shortly 
before the encounter and continues to drag material from the tidal 
arm (Figure [TJ even several kyr later. The slight decrease of mass 
of both structures is mainly caused by dynamical scattering and, 
in particular for the inner disc, accretion onto the star, but also 
by a small amount due the numerical diffusion that is unavoidable 
in SPH calculations. In a real disc the flnal masses would thus be 
somewhat larger. 



about 3000 yr. At the same time, the pre-existing disc suf- 
fers a radial compression (see frames b-d in Figure [1]). The 
distinction of these two structures has been done by the an- 
gular momentum grouping method described in Section [4. 11 
It has to be noted that the curves in Figure [J] show a repre- 
sentative fraction of each gas body defined by this method. 
Since part of the gas particles are also subject to scatter- 
ing into large separations (> 100 AU) and/or inclinations 
(> 15 degrees) these are omitted by this selection process, 
thus leading to a decreasing mass of each disc. These losses 
are in part due to dynamical scattering and frictional heat- 
ing of the gas, in part due to accretion by the bullet star. To 
a certain degree, they are also due to numerical diffusion, 
which is unavoidable in SPH, especially for relatively low 
resolution, and may thus be reduced by using larger particle 
numbers in future calculations. However, the mass-loss does 
not affect the central conclusion of this paper in any way. 



4.3 Pre-existing planets 

We also studied the effects of such a gas capture event on a 
pre-existing planetary system. Here, the circumstellar disc 
of the bullet star has been replaced by a model planetary 
system equivalent Jupiter, Saturn, Uranus, and Neptune. 
Like the stars, the planets are represented by sink particles 
with a sink radius of 0.5 AU. This value is a compromise 
between spatial resolution and mass resolution. Gas parti- 
cles passing through the sphere defined by the sink radius 
are considered as captured if they are gravitationally bound 
to the sink, otherwise they simply pass through it. The cur- 
rent results, summarised in Table [U show that a planetary 
system equivalent to the Sun orbited by Jupiter, Saturn, 
Uranus, and Neptune (with initial orbital radii of 5, 10, 20, 
and 30 AU, respectively) is severely affected by gas-planet 
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Table 1. IMasses and orbital elements of four model planets before 
and after a flyby and capture event equivalent to the scenario in 
Fig. [T] The planets, which replace the bullet star disc of the flrst 
scenario, are initially equivalent to Jupiter, Saturn, Uranus and 
Neptune. 20 kyr after the flyby, when the capture and accretion 
process is essentially over, the orbits are shrunk to a fraction of 
their initial values while being highly eccentric now. Furthermore, 
the orbital positions of planets 2-4 have been altered. The system 
has apparently become unstable. 



interaction. During the capture some 0.01 Mq of gas flow 
towards the bullet star, crossing the planetary orbits. In re- 
sponse, the planets migrate inward to semimajor axes of 
about one tenth their initial separations, i.e. between about 
0.5 and 4 AU. In addition, some orbits get strongly eccen- 
tric (up to e = 0.8, where e = means a circular orbit 
and e > 1.0 means ejection). Such a system is dynamically 
unstable, i.e. some planets are probably ejected over time, 
typically leaving one to three planets, with the most massive 
planet on a close-in orbit. Systems like these are expected 
to undergo significant subsequent evolution. Even after ejec- 
tion of the most unstable planets the remaining ones may 
still be subject to Kozai resonance, altering their eccentrici- 
ties and inclinations over many orbital times. But it is clear 
that these results suggest that misaligned planets and hot 
Jupiters as well as short-period planets on eccentric orbits 
may both be an outcome of gas flow onto a forming or ex- 
isting system. 



4.4 Stellar spin 

The current work assumes that the stellar spin is not largely 
changed since the assembly of its circumstellar disc, and thus 
its angular momentum vector points into approximately the 
same direction as the angular momentum vector of the disc 
before the encounte r. T his assumption, however, may not be 
true in general. As iLai et al. (2011!) have found, the stellar 
spin may, in some cases, be tilted through magnetic inter- 
action with the disc. Furthermore, during the early phases, 
when the protostar has still a considerable diameter (tenths 
of an AU), it may still re-align to the disc after an encounter 
event of the kind proposed in our paper. However, a large 
fraction of stars will already be compact enough to be ef- 
fectively decoupled from angular momentum transfer with 
their disc, given the absence of strong magnetic torque. Since 
encounter events and gas capture are a natural consequence 
of star formation in dense environments, we conclude that 
they indeed do play an important role in forming mis-aligned 
planetary systems. 



4.5 Resolution issues 

The required resolution strongly depends on the kind of 
scenario to be modelled. For example, for studying frag- 
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mentation the local Jeans mass has to be sufficiently re- 
solved, meaning by at least about 50 SPH particles. The 
study of the global behaviour of gas masses without focus 
on fragmentation, on the other hand, requires much less 
resolution. In our model, the mass of each gas particle is 
1.92-10"^ Mq or 0.64 Me. The well-resolved mass therefore 
is about 32 M® or twice the mass of Neptune. The trans- 
ferred mass in both models (encounter of disc with target 
and planetary-system vs. target, respectively) is of the or- 
der of 0.01 Mq or 5000 SPH particles, so mass resolution is 
not an issue here. For calculating the i mpact on dust com- 
ponents or to study vortex forma tion ([Barge fc Sommerial 
I1995I : iKlahr fc B odenheimerl l2004l l this resolution is likely 
to be insufficient. These issues will be addressed in future 
work. 

One could argue that the masses of Uranus and Nep- 
tune equivalents (planets 3 and 4, see Section [43)) . which are 
roughly equal to the interacting mass required for efficient 
momentum transfer, are less-optimally resolved by only 20 
to 25 SPH particles. The statistical scatter in terms of the 
number of gas particles interacting with these minor giant 
planets is about y^l/25 = 20 % and therefore does not sig- 
nificantly influence the outcome of drag-induced migration. 
Furthermore, this behaviour is fully consistent with the mi- 
gration of the Jupiter and Saturn equivalents (planets 1 and 
2). A word of caution has, however, to be stated at this point: 
Apart from mere resolution issues the proper treatment of 
accretion onto sinks shows a number of pitfalls. For example, 
the void of material in the sink may lead to artificial suction 
effects and thus speed-up accretion. The time- evolution of 
the planetary orbits has therefore to be treated with great 
caution. The long-term outcome, however, largely depends 
on basic physical principles like conservation of angular mo- 
mentum and can therefore be expected to reflect the real 
consequences on a pre-existing planetary system. 



5 DISCUSSION AND CONCLUSIONS 

We have analysed the general effects of inclined gas capture 
onto a pre-existing circumstellar disc upon passage through 
a dense gaseous reservoir. The secondary gas capture does 
not destroy a pre-existing disc even if the new material in- 
flows at a highly inclined angle. Rather, the inner regions 
of the disc (inside about 30 AU) become denser, while some 
of the outer disc material is scattered away. The major re- 
sult is the formation of an inclined combined circumstellar 
disc from captured material and the pre-existing disc. After 
capturing a few 0.01 Mq into an initially inclined annulus 
both structures tend to align with each other within a few 
10 000 years after the encounter. An annulus with an initial 
inclination of 135° aligns to about 35° within 20 000 years. 
Due to conservation of angular momentum this results in 
a net shift of more than 90° in our model, and even more 
for lower-mass pre-existing discs or disc-less stars. If planets 
form from the resulting inclined disc one naturally arrives 
at misaligned planets. 

The resulting disc may have two chemically distinct ra- 
dial regions, the inner region which is associated with the 
star, and the outer region composed from the captured ma- 
terial. 

In addition, the basic effect of capture onto a pre- 



existing planetary system have been estimated. Orbits typ- 
ically shrink and become strongly eccentric in response to 
the contact with captured gas and subsequent planet-planet 
interaction, while their orbital planes can be dramatically al- 
tered thus naturally leading to (even mutually) misaligned 
short-period planets with eccentric orbits. 

Important constraints are given by recent observations 
of apparent inclinations between circumstellar d iscs and 
their host stars. In a most recent analysis Wat son et al.l 
(|2011|) have found no significant misalignment between the 
normal vectors of debris discs and the spin axis of low-mass 
to solar- type stars. In particular, the projected inclinations 
between the star and its disc typically differ by 5-10 degrees, 
and 15-45 degrees in a few cases. Their results do not neces- 
sarily contradict our scenario since even moderate disc tilts 
and warps may lead to a Kozai resonance in a planetary 
system born out of such a disc. Fo l lowing the mechanism 
descri bed by Fabrvckv &: Tremaind (|2007[ ) and iNaoz et al.l 
(|201l|) near-random misalignment of a close-in planet and 
the spin of its host star may result. 

Future studies will also consider protoplanetary discs 
with embedded planets and systems with a brown dwarf 
on a wide orbit (about 100-200 AU). In addition, the conse- 
quences for dust content and mixing of differently composed 
material will be analysed. 

The mechanism described in this paper provides a nat- 
ural scenario for the formation of misaligned planetary sys- 
tems in gas-rich dense star-forming regions. While other 
models have been proposed our scenario requires rather sim- 
ple assumptions, and may even provide the initially mis- 
aligned planetary orbits required by the recent model of 
iNaoz et al. (2011). 
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